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Using the spectral energy distribution of M87, a nearby radio galaxy in the Virgo cluster, and assuming a 
supermassive black hole induced spike in the dark matter halo profile, we exclude any dark matter candidate 
with a velocity-independent (s-wave) annihilation cross-section of the order of (<jv) ~ 10~ 26 cm 3 s -1 and a 
mass up to 0(100) TeV. These limits supersede all previous constraints on thermal, s-wave, annihilating dark 
matter candidates by orders of magnitude, and rule out the entire canonical mass range. We remark in addition 
that, under the assumption of a spike, dark matter particles with a mass of a few TeV and an annihilation 
cross-section of ~ 10~ 27 cm 3 s -1 could explain the TeV y-ray emission observed in M87. A central dark 
matter spike is plausibly present around the supermassive black hole at the center of M87, for various, although 
not all, formation scenarios, and would have profound implications for our understanding of the dark matter 
microphysics. 

PACS numbers: 95.35.+d, 96.50.S-, 95.85.Pw 


I. INTRODUCTION 

Finding signatures of dark matter annihilations in the sky 
has always been a priority for the dark matter community. 
Not only would this validate the particle nature of dark mat¬ 
ter (DM) but it would also give some insights about dark 
matter properties. Evidence for an excess of cosmic-rays 
in a DM-rich environment (in particular our Galactic center) 
could provide indirect evidence for DM annihilations or de¬ 
cay. Evidence for a total velocity-independent annihilation 
cross-section of about (ov) ~3x 10“ 26 cm 3 s" * 1 would in ad¬ 
dition support the hypothesis that DM was once in thermal 
equilibrium with standard model particles. 

In the absence of annihilation signatures in DM halos, strin¬ 
gent limits are being placed on the DM self-annihilation cross- 
section as a function of the DM mass, wdm- The most severe 
limits originate from e.g. the diffuse y-ray background in the 
Milky Way and its companion galaxies (dwarf spheroidals), as 
well as from distortions of the cosmic microwave background 
(CMB) and both high energy positrons and antiprotons. Alto¬ 
gether these measurements already rule out the simples^] ther¬ 
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mal, velocity independent, dark matter scenarios with a mass 
ranging from a few MeV |2)[3] to ~ 100 GeV (see e.g. [4j- 
ITOl ). but here we show that one can go a step further. We use 
the spectral diffuse emission of M87, a nearby radio galaxy 
in the Virgo cluster located about 16 Mpc from us, to exclude 
heavier and more weakly interacting DM particles. 

We observe that the presence of a supermassive black hole 
(BH) in the core of M87 may increase the DM energy distribu¬ 
tion so much toward the galactic center that the predicted flux 
expected from thermal DM particles would exceed observa¬ 
tions. In the standard picture, the DM energy density follows 
a power law, p r 7 in the inner region, with r the distance 
from the galactic center and y ~ 1 the slope for a Navarro- 
Frenk-White (NFW) profile [TT]. However, very close to the 
black hole, the DM profile may rise very steeply. Assum¬ 
ing that the BH grew adiabatically, expectations are that y 
should instead lie between 2.25 and 2.5, with the typical value 
y = y sp = 7/3, according to Ref. II 1 21 . in the very inner region. 
Such an enhancement of the DM energy density is referred 
to as a spike and is expected to enhance the brightness of the 
electromagnetic flux originating from annihilating DM parti¬ 
cles by a factor (p S p/pNFw ) 2 ~ r 2 T~7sp) 5 a t least dose to the 
DM mass threshold. 

The existence of such a DM spike is however debatable. If 
the growth of the BH was not adiabatic (as expected if the BH 
seed were brought in by a merger), then the inner DM energy 
density profile would behave instead as p « r 4/3 Usui. 
Besides, if the DM halo itself underwent a merger, or if the 
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BH did not grow exactly at the center of the DM halo, the 
inner DM halo profile would follow p r ' ] > 2 , thus consid¬ 
erably reducing the electromagnetic flux expected from DM 
particles. Finally, even if a spike could form with y sp ~ 7/3, 
the process of dynamical relaxation by DM scattering off stars 
could smooth down the spike and lead to a DM halo profile of 
the form p °c r 3/ ’ 2 , which would correspond to a Moore pro¬ 
file fl5l . 

The latter argument may not be relevant for M87, though. 
Indeed dynamical heating by stars is inefficient when the dy¬ 
namical relaxation time t r in the core is larger than the Hubble 
time (~ 10 10 yr). This effect varies significantly from one 
galaxy to another, depending on the dynamical properties of 
the stellar core SBJ. M87 is dynamically young: its relaxation 
time is estimated to be t r ~ 10 5 Gyr (instead of ~ 2.5 Gyr for 
the Milky Way) due to the strongdependence on the velocity 
dispersion of the stars and DM. rl Hence, a spike formed at 
early times is much more likely to have survived galaxy dy¬ 
namics up to the present epoch in M87 than in the Milky Way. 

In what follows, we assume that a spike has formed in M87 
and, given the above argument, also survived the scattering 
off stars. We will study its impact on the electromagnetic 
signatures expected from DM annihilations and derive strin¬ 
gent limits on the DM properties. The paper is structured as 
follows. First we review the calculations of the electromag¬ 
netic flux from annihilating DM in Sec. [II] Then we derive in 
Sec.[III]the upper limits on the annihilation cross-section as a 
function of the DM mass in the presence or absence of a spike. 
In Sec. |IV[ we find the corresponding upper limits when one 
takes into account both a leptonic (or hadronic) jet and a dark 
matter spike in the inner part of M87. We further discuss the 
possibility of fitting the observed TeV y-rays with the prompt 
emission from DM annihilations in the spike. We conclude in 
Sec.[V] Technical details can be found in the Appendices. 

H. DIFFUSE EMISSION IN THE PRESENCE OF A DM 
SPIKE 

The spike is modelled by a DM energy density p °e r~ 7 / 3 , 
starting from the saturation radius, denoted by r sat and de¬ 
termined by the DM mass and annihilation cross-section (see 
Fig. [TV up to the spike radius r sp . Outside this inner region, 
i.e. for r > r sp , we assume a NFW profile ITO . More details 
are given in Appendix A. These assumptions are similar to 
those made in Ref. fT~8l —where we computed the boost in 
synchrotron radiation due to the presence of a spike in our 
galaxy—and to those of Ref. ESI where the authors studied 
the enhancement of the diffuse extragalactic y-ray background 
induced by DM spikes in other galaxies. In Fig. [T| we show 
the resulting profiles for two different values of the DM an¬ 
nihilation cross-section, so as to illustrate the impact of DM 
annihilations on the value of the saturation radius. As one can 


2 The relaxation time goes as t r o 3 mm, with G proportional to , 

where Mbh is the mass of the BH. 


see (and as is very well known), a larger cross-section tends 
to smooth down the profile. 
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FIG. 1. DM energy density as a function of the distance from the 
center, for a DM spike with y sp = 7/3. DM annihilations soften the 
central cusp differently, depending on the DM mass and cross-section 
(as illustrated by the dashed and solid lines for a fixed DM mass and 
two different values of the annihilation cross-section). 

The DM contribution to the spectral energy distribution 
(SED) of M87 is essentially three-fold. Photons can be pro¬ 
duced in the observed regime by (i) intermediate and final- 
state radiation from the charged particles which are exchanged 
or produced in the DM annihilation process, (ii) decay or 
hadronization of the particles in the final state, and (iii) inverse 
Compton scattering (ICS) of electrons produced by DM anni¬ 
hilations off low energy photons (CMB, infrared and UV), as 
well as synchrotron emission from these DM-induced elec¬ 
trons. We will neglect any possible bremsstrahlung emission 
since M87 is deficient in cold gas EOl . 

The first two emission modes (i and ii) are referred to as 
prompt emission. They are independent of the galaxy dynam¬ 
ics and are only determined by the DM particle physics prop¬ 
erties. On the contrary, the third type of emission—namely 
synchrotron and ICS—strongly depends on the properties of 
the interstellar medium in the galaxy, such as the magnetic 
field strength and interstellar radiation field. In the presence 
of a very strong magnetic field, synchrotron emission becomes 
the main source of low energy photons but it is also the main 
energy loss for the electrons and positrons produced by the 
DM. As a result, we find that prompt emission dominates ICS 
and is therefore the dominant source of high energy gamma- 
rays, while synchrotron radiation is the dominant source of 
X-rays. 

The assumption of a very strong magnetic field in this work 
is justified by the presence of a BH at the center of M87. We 
will assume typically B = 10 10 — 10 11 ,uG, corresponding to 
the equipartition model, as discussed in Refs. EQEI. As 
a result, the electrons and positrons that are produced by the 
DM annihilations in the inner region are expected to stay con- 
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fined to their site of injection, i.e. essentially in a sphere of 
radius r sp . This means that we can safely disregard spatial 
diffusion. 

We compute the prompt and synchrotron emission as fol¬ 
lows. The prompt diffuse y-ray intensity, / pl '° mp, (0), at angle 0 
from the center, can be estimated by integrating the DM halo 
density over the line of sight (l.o.s.) coordinate s: 


.prompt 


(9)=£y 2 


drij 

dEydQ. 


Ey (Cv),d Ny.j 

4jtr l m 2 DU d Ey 


I p 2 (r{s,B))ds 

J l.O.S. 


(1) 

where T| (= 2 here) is a factor that accounts for the nature 
(real/complex, Majorana/Dirac) of DM, dNyy/dEy is the y-ray 
spectrum that originates from the specific annihilation channel 
i (see Ref. (23]) and (or), the DM annihilation cross-section 
into this final state. Finally we use Ey = hv with h the Planck 
constant. 

To compute the synchrotron intensity, we need one ad¬ 
ditional step. We first need to determine the electron and 
positron spectrum from the DM annihilation rate, using 


Ve,i( r > £ ) 


1 (gv)i /p(r)\ 2 r »M d dE 
b{r,E ) r) V»dm/ Je d E s S 


( 2 ) 


and convolve it with the synchrotron power, P v (r,E) (see Ap¬ 
pendix B), so as to obtain the synchrotron emissivity 


rm DM 

jv,i(r)= 2 Pv{r,E)y eti {r,E)dE. (3) 

J m e 

The term dN e j/dEs in Eq. ([2]) represents the number of elec- 
tronf0produced by the decay or hadronization of the final state 
i. The total emissivity is the sum of the electron and positron 
contributions, hence the factor 2 in Eq. 0- to account for the 
fact that a positron is always produced simultaneously with an 
electron. 

The term b(r. E) corresponds to the total energy loss rate. 
Given the large values of the magnetic field that we consider in 
this paper, the losses are dominated by the synchrotron losses, 
that is 


b(r,E) =b syD (r,E) 


4 

3° TC 


B{r) 2 

2qo 


7 L , 


(4) 


where Ot is the Thomson cross-section, B(r) the intensity of 
the magnetic field, c the speed of light, yL = E/(m e c 2 ) the 
Lorentz factor of the electrons, m e the electron mass and qo 
the vacuum permeability, see e.g. Ref. |24l . The synchrotron 
intensity finally reads 



jv,i(r(s,Q)) 

4n 


ds. 


(5) 


Note that both injection spectra dNyj/dEy and dN e y/dE$ in 
Eqs. 0 and 0 are taken from Ref. l23l and include elec- 
troweak corrections as these become increasingly important 


3 d/Vy /d/y ~ 8(£s — »>dm) if the DM directly annihilates into electrons 
and positrons. 


at high energy. Also, to perform the line-of-sight integration, 
we will assume that the magnetic field distribution in the inner 
region is shaped by the accretion flow. Hence, unless stated 
otherwise, we will assume that the radial dependence of the 
magnetic field is given by the equipartition model and thus 
goes as B(r) r“ 5 / 4 , as discussed in Refs. [2l3 [22]. This 
specific form actually leads to very large values of the mag¬ 
netic field toward the center, typically up to 10 10 — 10 11 ,uG in 
the very inner region—as mentioned previously—which is at 
least eight orders of magnitude larger than the values usually 
considered in the Milky Way (see Appendix B). 

In principle, one should also take into account the effect 
of advection of electrons and positrons toward the center by 
the accretion flow around the BH, which increases the syn¬ 
chrotron flux typically in the range 10 12 — 10 14 Hz ]22l . How¬ 
ever, we disregard this effect throughout the present paper 
since most of our constraints come from higher frequencies, 
given the large magnetic field strengths we consider. More¬ 
over, even for smaller magnetic fields, including advection 
would not weaken our constraints but would only make them 
more stringent, so we remain conservative in this regard. 


III. UPPER LIMITS ON THE ANNIHILATION 
CROSS-SECTION 


Limits are set on the DM annihilation cross-section by com¬ 
paring the expected emission from DM with the measured 
SED for M87. Most data points have actually been compiled 
by the Fermi Collaboration in Ref. (25]. We use in particular: 

• the historical measurements of the core emission (from 
millimeter to X-rays (2614321 ). 

• the MOJAVE VLBA data point at 15 GHz which was 
derived in Ref. (251 (the data was reported in Ref. (33l ). 

• the 2009 X-ray data points which were derived in 
Ref. (251 from the 2009 Chandra measurements (34], 

• the 2009 Fermi-LAT data (25], 

• the 2004 HESS data (35], 

• the 2007 VERITAS data EQ, 

• the 2011 MAGIC data (37], 

which essentially give us the observed value of the electro¬ 
magnetic flux between 10 10 and 10 27 Hz. 


A. Methodology 

The dark matter contribution is estimated by integrating the 
prompt and synchrotron intensities v/ pr i ° mpl syn [see Eqs. (|T|) 
and 0] over a field of view that is centred on the galactic 
center and set by the angular resolution of the relevant exper¬ 
iment for a given frequency. Given the spherical symmetry of 










4 




FIG. 2. SED of M87 from the millimeter band to TeV y-rays, from a DM spike, for 100 GeV DM particles annihilating into with 
(ov) = 8.12 x 10~ 29 cm 3 s _1 (left panel), and into bb with (av) = 3.94 x 1CP 29 cm 3 s _1 (right panel). These are the maximal values of the 
cross-section compatible with the data for a DM candidate of 100 GeV. 


the spike, the prompt and synchrotron spectra are given by 

r 0 re s 

vF prompt,syn = 2jt ^ v /P™ mpt ’ Syn (0) S in0d0. (6) 

We recall that Chandra has an angular resolution of 0.5 arcsec 
over the whole energy range considered here l(38l while the 
angular resolution of the Fermi experiment reads 0.8° x E^y 
ll25l (ftGeV is the energy of the y-rays normalized to 1 GeV). 
HESS, VERITAS and MAGIC all have angular resolutions of 
the order of 0.1°, see Refs I39H4TI . In practice though, the 
spike is contained in such a small region that the exact value of 
the upper bound of the integral in Eq. |6) does not significantly 
affect the result. 

To set limits, we require that the synchrotron and prompt 
emission fluxes that are induced by the DM annihilations do 
not exceed the error bars on the flux for any measured data 
point. More specifically, we exclude any value of the annihi¬ 
lation cross-section that satisfies the following inequality for 
any observed frequency v: 

sr del - (S° bs + A S ° bs ) > K (S° bs + AS° bs ), (7) 

with S v = vf-\, and K < 1 (typically K = 10 4 ). The terms 
.S™’ del , S° bs and A S° bs represent respectively the expected DM 
contribution, the observed SED, and the la error bar at fre¬ 
quency v. 

Note that large values of the annihilation cross-section, i.e. 
(av) » 10“ 27 (»; D m/(10 GeV)) cm 3 s _1 , (8) 


flatten the inner part of the spike below a saturation radius 
given by 



( 9 ) 


Therefore one cannot always rescale the flux for dif¬ 
ferent values of (av), since in some cases the cross- 
section actually modifies the DM profile. For (av) < 
10- 27 (wdm/( 10 GeV)) cm 3 s _1 , on the other hand, the satu¬ 
ration radius is very small (below 10 2 pc) and in this regime 
the fluxes that we compute are simply proportional to the an¬ 
nihilation cross-section. 


B. Results 

In Fig. [2] we plot the largest allowed electromagnetic emis¬ 
sion (prompt plus synchrotron) expected from DM annihila¬ 
tions for a 100 GeV DM candidate. The left panel shows the 
predictions for DM annihilations into /j 1 /A while the right 
panel shows the predictions for annihilations into bb. The two 
bumps correspond to the synchrotron (left) and prompt (right) 
emission. We also derive the constraints on the annihilation 
cross-section for any DM mass and eight annihilation chan¬ 
nels (e + e~, T + TA qq, bb, tt, ZZ, hh. with li the standard 

model Higgs boson and q = u,d,s), shown in Fig. [3] The left 
panel shows the constraints in the presence of a spike and the 
right panel shows the constraints without a spike (assuming a 
NFW profile). 

The constraints in the case of a spike essentially rule out any 
DM candidate with a thermal (s-wave) cross-section, from a 
few GeV to a hundred TeVj^Jln other words, they rule out 


4 These constraints can be extended down to the MeV range for leptons and 

light quarks, provided the magnetic field is of the order of the equipar- 

tition value. Synchrotron emission then peaks around 10 1 °-10 11 * * * * GHz. 

For smaller values of the magnetic field, the synchrotron peak falls be¬ 

low 10 10 GHz, i.e. in the radio range, where synchrotron self-absorption 

significantly reduces the flux (see e.g. Ref. 1221 1. thus preventing one from 

setting any constraints on the DM annihilation cross-section. 
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FIG. 3. Upper limits on the annihilation cross-section as a function of DM mass for various annihilation channels. Constraints derived 
assuming a spike in the DM distribution are shown in the left panel while in the right panel we show the constraints obtained assuming a 
standard NFW cusp. 


the entire mass range relevant for thermal DM. The only ex¬ 
ceptions to our generic conclusions are for candidates which 
mostly annihilate into p 1 p (the limit is then about 30 TeV) 
or which annihilate democratically into all channels (but the 
limit would still be close to 100 TeV nevertheless). We also 
confirm that thermal candidates with a p-wave suppressed an¬ 
nihilation cross-section are ruled out if they are much lighter 
than a few GeV. For comparison, our limits in the case of a 
pure NFW profile are considerably weaker, see Fig. [3] (right 
panel). They only rule out very light (below a few GeV) s- 
wave thermal DM candidates. 

Of course our conclusions rely on the crucial assumption of 
the existence of a spike with y sp ~ 7/3. While this remains 
speculative, the stability of such a spike with respect to the 
scattering off stars is very likely. Hence, if the initial con¬ 
ditions were such that a spike could form in M87, our con¬ 
straints would rule out a very large chunk of the thermal DM 
parameter space. An alternative interpretation of our results is 
that the discovery of a thermal s-wave DM candidate would 
rule out the existence of a spike in M87. This would in turn 
constrain the evolution and formation of the supermassive BH 
at the center of M87. 


C. Robustness of our constraints 

In the previous sections, we have considered extremely 
large values of the magnetic field (several orders of magni¬ 
tude with respect to the Milky Way) and we have neglected 
absorption. It is therefore legitimate to question the robust¬ 
ness of our limits with respect to the magnetic field model and 
absorption processes. 


1. Dependence on the magnetic field 

Our most stringent limit on the DM contribution in M87 
is set by the Chandra X-ray data. Since most of the DM 
induced X-ray signal originates from synchrotron radiation 
and synchrotron emission strongly depends on the magnetic 
field, a weaker magnetic field could weaken our constraints. 
This is particularly worrying for the e' e and p p final 
states which give the largest X-ray contribution in M87 when 
m DM <100 GeV, as shown in Fig. [2] and which could become 
unconstrained. 

However a weaker magnetic field in the inner region, 
taken for example to be constant and about B = 10 5 pG> 
(as suggested in Ref. Il42l ). only weakens our constraints by 
about one order of magnitude below 30 GeV. We thus get 
(or) < 10- 29 cm 3 s _1 for e + e~ and p + p~ instead of (av) < 
10- 30 cm 3 s' 1 . Hence, even in the case of a weaker magnetic 
field in the inner region, we can rule out thermal s-wave DM. 
Note that decreasing the magnetic field to B = lmG makes our 
constraints stronger again as the signal would be constrained 
by the MOJAVE data. Finally, as noted in Ref. Il22l . if the 
magnetic field is significantly smaller than the equipartition 
value, synchrotron self-Compton emission decreases the DM- 
induced electron spectrum and thus also the synchrotron flux 
accordingly. Based on the results of Ref. l22l . we estimate 
that our limits for the e' e and p 1 p channels are weakened 
by an additional order of magnitude below ~ 100 GeV for 
a magnetic field strength weaker than 10 5 ,uG. However this 
does not affect our conclusion since we can still exclude ther¬ 
mal s-wave DM. 

A magnetic field of the order of 10 5 pG also changes the 
limits for candidates with a mass above ~ 50 TeV for the e + e~ 
and p p channels. Indeed, the synchrotron peak then falls in 
the energy range [100 keV, 100 MeV] where there are no data. 
In that case, the limit is given by the prompt component and 
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weakened to the level of 3 x 10“ 25 cm 3 s'. As a reminder, 
in the case of a stronger magnetic field, we could rule out the 
canonical thermal cross-section for the e + e~ channel. Thus, 
for such a value of the magnetic field, one can no longer ex¬ 
clude s-wave DM for particles heavier than 50 Te V. It is worth 
pointing out though that a smaller value of the magnetic field 
of e.g. 10 3 ,uG would not alleviate our constraints as it would 
give rise to an excess in X-rays, which has not been observed. 
Such a value would therefore lead to an exclusion limit in¬ 
stead, similar to the one obtained for the equipartition mag¬ 
netic field. 

Finally, the annihilation channels which give a softer elec¬ 
tron spectrum, e.g. the bb channel, are unaffected by a weaker 
magnetic field since the prompt y-ray emission dominates the 
synchrotron emission, see Fig. [2] and therefore our limits for 
these channels are independent of the magnetic field. 
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2. Absorption 

Absorption is another process which could weaken our con¬ 
clusions. We may have overestimated the flux by not account¬ 
ing for the photons which have been emitted by synchrotron 
radiation and absorbed by the same electron population that 
produced them. Since the authors of EU [22l showed that 
this effect is only very efficient below 10 1() Hz, we cut the 
synchrotron emission below this critical frequency. This pre¬ 
vents us from constraining the scenarios emitting in this en¬ 
ergy range, i.e. candidates lighter than 0{ 1) GeV, unless the 
magnetic field is strong. 

At the other end of the spectrum, high-energy y-rays could 
also be absorbed via e + e~ pair production with the back¬ 
ground radiation field. However, the authors of Refs. (421331 
showed that the inner region of M87 is transparent to y-rays 
below ~ 10 TeV. Since we only have data below 10 TeV and 
absorption is relevant only above 10 TeV, we can neglect this 
effect for all candidates below 10 TeV. For the much heavier 
candidates, the lack of data above 10 TeV makes absorption 
irrelevant for the moment. Hence we have neglected absorp¬ 
tion in our study. 

IV. DM SPIKE AND JET 

In the previous section, we have investigated the DM con¬ 
tribution to the SED of M87 but neglected the contribution 
from the BH. In reality the jet emission associated with the 
BH must be taken into account. Indeed, to be observable, any 
putative emission from DM should be brighter than the emis¬ 
sion from the jet. 

A. Jet emission 

The mechanism giving rise to high energy photons from the 
jet is not well known. It is unclear whether these photons have 
a leptonic or hadronic origin. In the most popular model, the 
y-rays originate from electrons contained in a blob of plasma 


FIG. 4. SED of M87 from the millimeter band to TeV y-rays. The 
orange solid line represents the SSC model that best fits the data, for 
a 55 mG magnetic field in the plasma blob, a Doppler factor of 3.9, 
and a radius of 4.5 mpc for the blob in its rest frame, as found in 
Ref. (25]. The blue dashed line corresponds to the signal expected 
from DM-jet scattering, as described in Ref. jH. Details on the 
parameters can be found in the text. 


that moves relativistic ally and possesses a fairly strong mag¬ 
netic field. This leptonic synchrotron self-Compton emission 
(SSC) proceeds in two steps: the electrons from the jet first 
produce photons in the infrared band due to synchrotron radi¬ 
ation in a strong magnetic field. Then, in a second step, these 
low energy photons are upscattered to y-ray energies by ICS 
on the same electron population that produced them. 

To take this effect into account and fit the spectral energy 
distribution of M87, we use the model described in Ref. |45) 
and take the best-fit SSC parameters given in Ref. |25l . see 
Appendix C. The parameters we consider are: a Doppler fac¬ 
tor So = 3.9, a magnetic field B = 55 mG, and a source radius 
R' b =4.5 mpc in the rest frame of the blob. The data that 
have been used for the fit are the 2009 MOJAVE, Chandra and 
Fermi-LAT data. The best-fit value for the normalisation of 
the electron distribution is K = 5.81 x 10 51 . The correspond¬ 
ing SSC emission for this set of parameters is shown in Fig. [4] 
(see the orange solid line). 

Note that the scattering of the DM particles off electrons 
and protons in the jet might also produce high energy photons 
(44ll and, consequently, lead to a characteristic signature in 
the Fermi-LAT data. The associated flux is proportional to the 
integral of the DM density over the line of sight (8 dm. as in 
B31), and the jet power (L). For M87, the highest allowed jet 
power is L ~ 10 45 erg s 1 l25l . For an optimal configuration 
of the DM spike, i.e. for the largest possible DM energy den¬ 
sity (correspondingly to the smallest possible saturation ra¬ 
dius, typically (cv) ~ 10~ 30 — 10~ 29 cm 3 s” 1 for m dm ~ 100 
GeV), the line-of-sight integration gives 8 dm ~ 10 9 M 0 pc 2 
which leads to a y-ray flux of ~ 10~ 14 erg cm~ 2 s . This 
is roughly two orders of magnitude below the Fermi data, as 













7 


shown in Fig. [4] (blue dashed line). This process is therefore 
subdominant for M87, and we will disregard it in the follow¬ 
ing discussion. 



FIG. 5. Upper limits on the annihilation cross-section as a function 
of DM mass for usual annihilation channels, obtained after summing 
the DM contribution to the photon emission expected from the jet 
(using the SSC model) and excluding cross-sections that depart from 
the best fit at 2 a. 


B. Upper limits on the annihilation cross-section with 
spike+jet 

Because the jet emission associated with the SSC model 
fits the data very well up to y-ray energies of 100 GeV, there is 
little room for a dark matter contribution to the SED of M87 
for candidates lighter than 20 TeV. Our results are shown in 
Fig-0 We have excluded any values of the annihilation cross- 
section that cause the total flux to depart from the best fit to 
the Chandra and Fermi data by more than 2c and improved 
the limits by about one order of magnitude with respect to 
the constraints derived without considering the jet emission. 
Above > 100 GeV, a DM contribution improves the fit and 
prevents us from setting a better limit. Note that we can also 
exclude thermal p-wave DM up to ~ 40 GeV. 

Because the jet emission fits the whole spectrum up to 100 
GeV and, in particular, fits the Chandra data which have the 
smallest error bars, any additional DM contribution—even 
small—tends to worsen the chi-square statistic and thus leads 
to stronger constraints. Since the Chandra data constrain the 
synchrotron contribution which is very sensitive to the mag¬ 
netic field, these constraints strongly depend on the strength 
of the magnetic field. 

Finally, let us recall that the limits derived in this section 
depend on the underlying jet model which, as we mentioned, 
is still very uncertain. Yet, the limits of Fig. [5]do illustrate the 
importance of including a model for the jet. 


C. Explaining the TeV data with a DM spike 

1. Fits with a jet + DM spike 

As shown in Fig. 0 the simplest SSC model does not 
explain the TeV emission measured by HESS, MAGIC and 
VERITAS (although some of the points are consistent with 
the jet model)0This led the authors of Ref. Ii47l to discuss the 
possibility that prompt emission from TeV DM may alleviate 
the discrepancy between the jet model and the TeV data. As¬ 
suming a NFW profile and the presence of DM clumps, they 
fit the data with a very large value of the annihilation cross- 
section (typically 3 x ICC 24 cm 3 s -1 ) and a very large boost 
factor of almost 1000. However, as shown in Fig. 0 in the 
presence of a spike we can fit the TeV data for a value of the 
annihilation cross-section smaller than the thermal value and 
no additional boost factor is needed. 

Using the Chandra, Fermi-LAT, VERITAS, MAGIC and 
HESS data (namely 24 data points), our best fit for the jet 
emission model gives % 2 « 85.5, that is % 2 / d.o.f. ss 3.56 for 
24 degrees of freedom. Adding a DM contribution (mod¬ 
eled by 2 free parameters, namely the DM mass and anni¬ 
hilation cross-section) to this best fit background model con¬ 
siderably improves the quality of the fit. We obtain yj ss 29.4 
(X 2 /d.o.f. « 1.34 for 24 — 2 d.o.f) for the bb channel (Fig. 0 
left panel) andx 2 «28.7 (x 2 /d.o.f. « 1.30) forthei+ir chan¬ 
nel (Fig.0 right panel). The corresponding best-fit values for 
the mass and cross-section are given in Table [T] 

The associated la and 2a confidence contour plots for both 
channels are shown in Fig. 0 For completeness, we also in¬ 
dicate the best-fit values obtained for other annihilation chan¬ 
nels. Note that the cc and gg channels are degenerate with the 
qq channel, so the same conclusions apply. Similarly W W 
and ZZ are also degenerate. 

2. Dependence on the magnetic field 

The best-fit values displayed in Table [I] are obtained by as¬ 
suming the same magnetic field intensity as in Sec. [II] and 


TABLE I. Best-fit DM mass and annihilation cross-section, for vari¬ 
ous characteristic annihilation channels. 


Channel 

m dm (TeV) 

(ov) (cm 

V 

- 1 ) 

X 2 /d.o.f. 

bb 

23+g 6 

3.9tf% x 

10 

-27 

1.34 

rl 

+ 

c-| 

1 

2 l +0 - 7 
z -0.5 

4 4+ 1 - 2 x 
—0.9 x 

10 

-28 

1.30 

qq 

16 + 1 4 

2.7±tl x 

10 - 

-27 

1.46 

It 

31+ 24 

6.0 + 4 3 x 

10 - 

-27 

1.33 

ZZ 

18 + 1 4 

3 - 9 -L5 x 

10 - 

-27 

1.29 

hh 

22^ 5 

4.3± 2 j x 

10 - 

-27 

1.25 


5 Also according to the authors of Ref. ®, it may be possible to refine the 

SSC model so as to fit the TeV data. 
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v (Hz) v (Hz) 


FIG. 6. SED of M87 from the millimeter band to TeV y-rays. The SSC model for the jet gives a double peak structure (orange solid line). The 
contribution from the DM spike is depicted by the purple dot-dashed line, for annihilations into bb (left panel) and t + t~ (right panel), with 
the synchrotron peak around 10 22 Hz and the prompt emission peak around 10 26 Hz. The black dashed line is the total SED. 



m DM ( Te V) 

FIG. 7. Confidence contours at the lo and 2o levels, in the plane 
annihilation cross-section vs DM mass, for the bb channel (black 
thick contours) and the channel (orange thin contours). The 

best fit points at the center of the contours correspond to ojdm = 
23 TeV, (ctv) = 3.9 x 1(T 27 cm 3 s _1 for bb, and m DM = 2.1 TeV, 
(ov) =4.4x 1(T 28 cm 3 s _1 fort+x - . 


reflect the fact that very heavy DM candidates give rise to 
very high energy y-rays. However this statement depends 
on the magnetic field and heavy particles can emit light at 
much lower energies. For example, we observe that a DM 
candidate with wdm ~ 20 TeV can lead to an excess of X- 
rays if the magnetic field is relatively weak (typically about 
10 5 — 10 6 /rG) and be ruled out by the Chandra data. 

Whether a candidate is ruled out or not, however, also de¬ 
pends on the model for the jet emission. By varying both the 
jet model and the DM component we can, for example, rec¬ 
oncile a DM candidate with ;«dm ~ 20 TeV (supposedly ruled 


out by the Chandra data in the presence of a relatively small 
magnetic field) with a possible noticeable contribution at TeV 
energies. Note that for such relatively small values of the mag¬ 
netic field, ICS and SSC are still negligible. 

If the magnetic field is even smaller, typically ~ 10 3 ,uG in 
the inner region, the synchrotron emission gives a signature 
at energies of a few eV corresponding to frequencies of about 
~ 10 15 Hz. In that case there is no tension with the X-ray data. 
However, ICS becomes non-negligible for moderate magnetic 
fields. The SSC emission could also be important but we ex¬ 
pect it to be subdominant. For B ~ 10 3 ,uG and wdm ~ 20 TeV, 
we expect ICS to give an additional contribution at TeV en¬ 
ergies, thus strengthening the case for an explanation of the 
observed high energy emission in terms of DM. 

All these remarks show that the best-fit values obtained by 
fitting prompt y-ray emission give a very good estimate of 
the contribution of annihilations from a DM spike to the TeV 
emission, fairly independently of the magnetic field and in¬ 
terstellar radiation field model. Therefore, our conclusion is 
that if there is indeed a DM spike in M87, then the subsequent 
annihilations can account for the TeV y-ray emission, with 
annihilation cross-sections 10 times smaller than the thermal 
value or even smaller depending on the channel. 


V. CONCLUSION 

We believe that the case for a DM spike at the center of the 
M87 galaxy is very strong. One therefore expects a signifi¬ 
cant annihilation signal from thermal DM candidates. In this 
paper, we have confronted the observed SED of M87 with the 
predicted emission from DM, and set extremely strong upper 
limits on the annihilation cross-section of DM particles as a 
function of the DM mass. These limits exclude thermal DM 
candidates with a velocity-independent (s-wave) cross-section 
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and a mass up to 0(100) TeV. Our results are independent of 
the magnetic field distribution and absorption processes what¬ 
ever the DM mass; the sole exception is for annihilations into 
light leptons but our conclusion remains valid for DM masses 
up to 50 TeV. Also we have shown that in the presence of 
a DM spike, TeV DM can explain the TeV y-ray data for 
annihilation cross-sections smaller than the canonical value 
(~ 10- 27 cm 3 s" 1 ). 

We expect similar constraints for galaxies containing a su- 
permassive BH with the same mass as in M87. Should such 
a spike be found, for example using stellar kinematics, one 
would exclude a very large chunk of the thermal DM param¬ 
eter space. This opens up a new path in DM searches, with 
great potential to elucidate the nature of DM particles. These 
results provide a strong motivation to look for further evidence 
for DM spikes in galaxies. 
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APPENDIX A: DM SPIKE MODEL 


rlO 5 

J4R S 


10 5 /? s 


47tr 2 p(r) dr = AM B h, 


(13) 


with AMbh = 5 x 10 8 M 0 . Considering the complex de¬ 
pendence of p on po, we use the fact that the mass is 
dominated by the contribution from r R$, i.e., typi¬ 
cally r > R m [ n = 0( 1 GOA'S). In this regime we have p ~ 
p S p(r). We can also factorize the dependence on po in 

Ps P , p S p(r) = syMPo~ T (K P / r oy y {K v / r V SV ’ with K P = 
a y ro (MbhAo) 3 Y ’ anc * we finslly obtain 


Po = 


4kR, 


(3 — y S p) AMbh 

As P -y r y ( £>3—Y S p _ /? 3-y sp 


4—y 


(14) 


sp 


'0 \ R 0 


We consider the profile derived by the authors of Ref. El 
for a DM spike growing from an initial profile °< po (r/ro) 7 : 


0 


p(r) = < 


Psp(r )Psat 


Psp(r) T Psat 



-2 


r < 4J?s 

4Rs <r<R sp 

r > R sP , 


GO) 


where the saturation density determined by DM annihilations 
reads 


Psat — 


mpM 
(av) /bh ’ 


( 11 ) 


with mp> m and (av) respectively the mass and annihilation 
cross-section of the DM particle, and ?bh the age of the BH. 
The spike profile reads 


R \ Isp 

Ix sp ' 


( 12 ) 


Numerically, we get po « 2.5 GeV cm 3 for y= 1. Finally, 
the total mass within 50 kpc is ~ 4 x 10 12 M e , compatible 
with the value derived from observations, 6 x 10 12 M e 1®. 

For completeness, we also consider the case of a DM cusp 
without a spike. In that case the profile is given by: 


PM 


[° 

J Psat 



r < 4R$ 

4Rs ^ f < Kat (15) 
r > r s at, 


where r sat = ro (po/p sa t) v , with the same value of po as in the 
presence of a spike. 

In practice throughout the paper we take y = 1, which cor¬ 
responds to the NFW profile ifTTIl . The corresponding spike 
has a power-law index of y sp = 7/3. 


Psp(r) = p R g y (r) 


r 






10 


APPENDIX B: SYNCHROTRON AND PROMPT EMISSION 
INTENSITIES 


Unless otherwise stated, we used the equipartition magnetic 
field model, characterized by: 


B(r) = 



r < r acc 

r a cc < r < r c 
r > r c . 


(16) 


We take Bq = 10 ,uG for the large scale value of the mag¬ 
netic field outside the inner cocoon of radius r c ~ 10 kpc seen 
for instance by LOFAR |49| . The authors of Ref. l2ll esti¬ 
mate the radius of the accretion region as r acc = 2 GMbh/va 0W , 
where Vfl ow ~ 500 — 700 km s is the velocity of the Galactic 
wind at the center of the Milky Way. Here we assume sim¬ 
ilar characteristics for the wind at the center of M87, so we 
just rescale the BH mass. For the Milky Way, the size of the 
accretion region was ~ 0.04 pc. Now, considering that the 
black hole in M87 has a mass approximately 1.5 x 10 3 times 
larger than the one at the center of the Milky Way, Sgr A*, we 
estimate r acc ~ 60 pc. The resulting equipartition magnetic 
field can reach very large values at the center, typically up to 
10 10 — 10 11 ,uG in the very inner region. 

To compute the synchrotron intensity, we first need the 
electron (and positron) spectrum, given by (see, e.g., ED) 


Ve,i( r > £ ) = 


1 (cv),- ( p(r) 


b(r,E) q \' w dm/ Je dE s 


f'”DM dN c 


(E S )dE s , 


(17) 

where we use q = 2 —which corresponds to the assumption 
that the DM particle is a Majorana fermion—and d/V e .,/tl/:s is 
the electron or positron injection spectrum that we take from 
Ref. l23l for each channel denoted by i. We use injection 
spectra that include electroweak corrections that become very 
important for large masses. b(r,E) is the total energy loss rate. 
Since we consider large magnetic fields in the inner region, 
essentially all the energy of electrons and positrons is lost in 
the form of synchrotron radiation and the contributions from 
ICS and synchrotron self-Compton turn out to be negligible 
ll22l . so that the loss term reads (see, e.g. El) 


b(r,E) = b sya (r,E) = (18) 

3 2q 0 

where Ox is the Thomson cross-section, B(r) the intensity of 
the magnetic field, c the speed of light, Yl = E/ ( m c c 2 ) the 
Lorentz factor of the electrons, m e the electron mass and ,uq 
the vacuum permeability. 

Then the synchrotron emissivity reads (see e.g. Ref. ED) 
r m dm 

jvA r ) = 2 PvM)Ve,M)d£, (19) 
Jm e 


chrotron emission spectrum is given by (see, e.g.. Ref. El) 


P,(r,E) 

47t£o W e C 


V 

Vcpj 


( 20 ) 


where e the elementary charge, £0 the vacuum permittivity, 
and the critical frequency is given by 


, N 3 eE 2 B{r) 
Vc ^ £ )- 4K,nh A ' 


( 21 ) 


G; is the isotropic synchrotron spectrum, obtained by averag¬ 
ing the synchrotron spectrum over an isotropic distribution of 
pitch angles 1(241 : 

Gi(x) = -[ g(— —)sin 2 ada, (22) 

' 2 Jo Vsina/ 

with G(f) = t f t °°Ks /3 (u) d u, where K$ / 3 is the modified Bessel 
function of order 5/3. To simplify the numerical treatment of 
the angle average, one may use, for instance, the parametriza- 
tion described in Ref. Bill . 

From there, the specific intensity for a given angle 0 from 
the center is given by the integral of the emissivity over the 
line of sight (l.o.s): 


rsyn 

'v,i 



./v,,(r(s,9)) 

47t 


d,v. 


(23) 


s being the radial coordinate along the line of sight and 
r(s, 0 ) = \Jd 2 + s 2 — 2ds cos Q & yj{d — s) 2 +dsQ 2 . The ap¬ 
proximation of small angles is justified since the characteristic 
radius Rm 87 of M87 (typically 50 kpc) is much smaller than 
the distance of M87, d = 16 Mpc. Also in practice we perform 
the integral over the l.o.s. between d — R\w an d d + Rm87i 
considering the concentrated nature of the DM profile. 

The specific intensity for prompt y-rays is simply given by 
the integral over the line of sight of the DM density squared 
(see, e.g.. Ref. |52| and references therein), 


.prompt 

VI v,i 


(Q)=E 2 


d nj 

d£j,dD 


E 2 (av)f d Nyj 

47 tq ml M dE y 


[ p 2 (r(s,0)) ds, 

Jl.o.s. 

(24) 


where dNyj/dEy is the prompt gamma ray spectrum taken 
from Ref. Il23ll and Ey = /;v with h the Planck constant. 
The specific intensity at 1 TeV as a function of the angle 
from the center, for a DM spike with y sp = 7 /3, for anni¬ 
hilations proceeding to bb, with jhdm_= 23 TeV and (ctv) = 
3.9 x 10- 27 cm 3 s 1 , is shown in Fig. ml The integral over the 
l.o.s. is approximately constant at the center, due to the van¬ 
ishing density below 4Rs, hence the plateau below ~ 100 /./as. 
Above ~ 1 arcsec, the change in slope is related to the outer 
part of the DM profile, assumed to follow the NFW distribu¬ 
tion. 


APPENDIX C: SYNCHROTRON SELF-COMPTON MODEL 
FOR THE JET 


where the factor 2 refers to the fact that both an electron and Here, we summarize the leptonic SSC model for the spec- 
a positron are produced in one DM annihilation, and the syn- tral energy distribution of M87 described in Ref. ED and 
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r (pc) 

10' 3 10' 2 10' 1 10 ° 10 1 10 2 10 3 10 4 



where y = 4 x 10 3 is the Lorentz factor at the break and Y> = 
10 7 is the maximum Lorentz factor of the electrons. 

From there, the observed SSC flux is given by m 


fSSC 
j Es 


9aTe' s 2 





(28) 


where f v , min = R(,/(8 dc) ' s the variability time scale of the 
source, R' b being the (comoving) radius of the blob, e' is the 
dimensionless energy of the scattered photon. The ICS pro¬ 
cess is encoded in Fc(g,r) which reads 


F c (q,r)=2qlnq+(l+2q)(l-q)+ ^ (1 -?) (29) 

if 1/(4Y 2 ) < q < 1 and Fc{q,T) = 0 otherwise, q and T are 
given by: 


FIG. 8. y-ray specific intensity at 1 TeV as a function of the angle 
from the center, for a DM spike with y sp = 7/3, for the bb channel, 
m DM = 23 TeV and (av) = 3.9 x 10~ 27 cm 3 s -1 . 


used by the Fermi Collaboration in Ref. {25l . All the primed 
quantities are defined in the rest frame of the plasma blob. 
Considering that the redshift of M87 is 0.00428£]we neglect 
redshift effects in our discussion. Following the notations 
of Ref. l45l . the observed synchrotron flux vF v is denoted 
/ E , where 8 = hv/(m e c 2 ) is the dimensionless energy of the 
emitted synchrotron in the observer’s frame. Similarly to the 
synchrotron flux obtained in Appendix B, the observed flux, 
rewritten in terms of dimensionless energies, and taking into 
account the Doppler boost, reads: 


q = 


e's/Y 


r ( l - e '/ Y )’ 


r = 4e'Y. 


(30) 


The kinematically allowed range of values for q translates into 
the integration bounds in Eq. (|28|: 



(31) 


Ym 


e' - e' 

y 2 


e' > e' 
e' < e'. 


(32) 


syn __ \/35Ye'e 3 Bc 
4jt hd 2 


j x A 7 e(Y)G i(x(e',Y)), 


(25) 


where 8 d is the Doppler factor, e' = e/So, and d = 16 Mpc is 
the distance of M87. N r e = n' e V b is the electron distribution in 
the rest frame of the blob, with n' e the electron number density, 
and V b ' the volume of the blob, x = v'/v' is rewritten in terms 
of the dimensionless quantities: 


x(e',Y) 


47te'/n 2 c 2 
3 eBhy 2 


(26) 


Note that the intensity of the magnetic field B is not primed but 
is also defined in the rest frame of the blob. For 6’,, we take 
the parametrization of Ref. ED, also used in Ref. H31 . For 
the electron distribution, we consider as in Ref. f25j a broken 
power law: 


N'M)=K 


y-pi 

yP2-piy-p 2 


i < Y < Y 

Yi <Y<Y 2 , 


( 27 ) 


6 http://messier.seds.org/xtra/supp/m_NED.html 
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